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Abstract 

Solvation  states  of  single  solvent  electrolyte  systems  of  ethylene  carbonate  (EC),  propylene  carbonate  (PC),  dimethyl  carbonate  (DMC), 
ehylmethyl  carbonate  (EMC )  and  diethyl  carbonate  (DEC )  with  LiPF6  were  characterized  by  ,3C-NMR  solvation  shift  and  molecular  orbital 
(MO)  simulation.  Dissociation  constants  and  solvation  constants  were  estimated  by  parameter  fitting  to  solvation  shift  using  a  simple 
equilibrium  model.  The  solvation  shifts  A  8  were  observed  not  only  at  a  lower  field  but  also  at  a  higher  field  due  to  change  of  net  charge  A p 
in  solvent  molecules  by  Li+  attachment.  This  particular  feature  of  solvation  shifts  was  demonstrated  in  the  molecular  orbital  simulation  as 
driven  by  the  change  of  net  charge  using  a  1:1  (Li + : solvent)  solvation  model.  ©  1997  Elsevier  Science  S.A. 
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1.  Introduction 

The  organic  solvent  chosen  determines  the  physical  prop¬ 
erties  of  the  electrolyte  and  cell  performance  of  a  lithium-ion 
secondary  battery.  However,  the  choice  of  the  electrolyte  is 
still  a  trial-and-error  process  because  even  solvents  having  a 
similar  chemical  structure  show  a  different  behavior.  To 
design  an  electrolyte  by  molecular  simulation,  it  is  necessary 
to  correlate  the  physical  properties  of  the  electrolyte  with  the 
chemical  structures  of  solvent  and  anion.  Recently,  molecular 
orbital  (  MO)  simulations  have  been  applied  to  elucidate  the 
solvation  structure  [  1  ] .  This  work  aims  at  understanding  the 
solvation  states  of  cyclic  and  linear  carbonate  compounds  in 
a  single  solvent  electrolyte  system. 

2.  Experimental 

Commercially  available  organic  solvents,  ethylene  carbon¬ 
ate  (EC),  propylene  carbonate  (PC),  dimethyl  carbonate 
(DMC),  ethylmethyl  carbonate  (EMC)  and  diethyl  carbon¬ 
ate  (DEC)  and  LiPF6  were  used.  The  0.5,  1.0  and  2.0  M 
single  solvent  electrolytes  were  prepared  in  an  argon-filled 
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glove  box.  ,  3C-NMR  chemical  shifts  of  organic  solvents  and 
electrolyte  were  recorded  on  a  JEOL  JNM-GX400  FT-NMR 
with  an  external  standard  of  tetramethylsilane  (TMS)  at 
25  °C,  except  for  EC  which  was  recorded  at  40  °C. 

Solvation  shifts  (A 8)  were  estimated  by  Eq.  (1),  where 
<5S  and  8{}  represent  chemical  shifts  of  electrolyte  (with  salt) 
and  solvent  ( without  salt),  respectively 

A  8=8s-80  (1) 

A  8  is  a  function  of  the  molar  ratio  (r)  of  LiPF6  to  the 
solvent  and  is  denoted  A<5(r).  Because  the  chemical  shifts 
of  electrolyte  are  observed  as  an  average  of  free  solvent  and 
solvated  solvent,  the  critical  solvation  shift  A  8'  ( the  net  shift 
for  solvated  molecule )  must  be  estimated  by  analyzing  the 
experimental  A  <5(  r).  Using  the  simplest  equilibrium  process 
written  in  Eqs.  (2)  and  (3),  the  equilibrium  constants  of  Kd 
(dissociation  constant,  inverse  of  association  constant  Kd)% 
( solvation  constant)  and  A  8 '  were  estimated  by  parameter 
fitting  to  the  experimental  A  8(  r)  [2] 

Ad 

LiPF6  Li +  +  PF6“  ( 2 ) 
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Li +  +  solvent  Li +  (solvent)  (3) 


J.  Aral  et  a l.  /  Journal  of  Power  Sources  68  ( 1997)  304-306 


305 


3.  Computational  method 

Molecular  orbital  calculations  were  performed  with  the 
MOPAC  version  94. 1  program  on  an  IBM  workstation  asso¬ 
ciated  with  CAChe^1  software.  A  1 : 1  ( Li +  :solvent)  solvation 
structure  model  binding  with  carbonyl  oxygen  was  optimized 
with  and  without  the  COSMO  {conductor-like  screening 
model)  method  [3].  Change  of  net  charges  on  the  solvent 
molecule  ( Ap)  was  estimated  by  comparing  the  calculation 
results  for  a  solvent  molecule  in  the  neutral  state  (p0)  and  the 
solvation  state  (ps),  as  expressed  by  Eq.  (4).  The  solvation 
enthalpy  A  H  was  calculated  by  Eq.  ( 5 ) ,  where  Hs  is  the  heat 
of  formation  for  the  solvation  state,  H0  is  that  for  the  neutral 
state  and  Hu  is  that  for  Li + . 

Ap  =  ps-p0 

\H  =  H,-(H0  +  HLi)  (5) 

4.  Results  and  discussion 

Fig.  1  shows  the  solvation  shift  A  §  as  a  function  of  the 
molar  ratio,  r,  along  with  the  calculated  curve  determined  by 
parameter  fitting  for  the  PC  electrolyte  system.  As  shown  in 
this  figure,  the  solvation  to  Li+  caused  not  only  a  lower  field 
shift  (charge  reduction),  but  also  a  higher  field  shift. 

Table  1  summarizes  Kd ,  Ks  and  critical  solvation  shift  for 
A  8'  each  carbon  atom  of  the  solvent  molecules.  Higher  values 
of  Kd  and  Ks  for  EC  may  be  due  to  the  higher  observation 
temperature  ( 40  °C ) .  Kd  and  ATS  for  EMC  were  highest  among 
linear  carbonate  compounds,  suggesting  a  steric  effect  due  to 
an  asymmetric  structure.  Kd  and  K s  for  DEC  were  lower,  this 
must  be  due  to  its  large  alkyl  group.  >  C=0  had  a  larger  A  8f 
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Fig.  1.  Solvation  shift  as  a  function  of  salt  to  solvent  molar  ratio  for  PC 
electrolyte  system.  The  calculated  curves  are  drawn  with  Kd,  and  the 
critical  shift  determined  by  parameter  fitting. 

in  the  cyclic  carbonate  compounds,  while  CH2  carbon  showed 
a  larger  A  8f  in  the  linear  carbonate  compounds.  A  8'  of  carbon 
beside  -O-  must  be  large  in  the  case  of  solvation  using  -O- 
oxygen.  This  suggests  that  >C=0  oxygen  bonding  was 
dominant  in  cyclic  carbonate  compounds,  while  -O-  oxygen 
bonding  may  take  place  in  linear  carbonate  ones. 

Table  2  lists  enthalpies  A H(e)  calculated  with  the 
COSMO  method  using  the  dielectric  constant  of  the  pure 
solvent,  enthalpies  calculated  in  vacuum  A//(  1 )  and  Ap(e) 
and  Ap(  1)  using  the  >C=0  bonding  1:1  solvation  model. 
Ap(  1 )  could  explain  the  sign  of  A  8'  but  not  in  good  order. 
Thus,  the  sign  change  of  A  8 '  was  due  to  charge  alternation 
of  solvent  molecule  caused  by  bonding  to  Li +  as  illustrated 
in  Fig.  2.  Ap(  1 )  for  >  C=0  was  much  larger  than  those  for 
CH3  and  CH2.  While  A  p(e)  did  not  give  negative  values  for 
the  CH3  group  in  PC  and  DEC,  it  gave  no  positive  values  for 
>  C=0  carbon  in  PC  and  EC.  There  are  two  ways  to  explain 
these  results:  (i)  e  of  the  electrolyte  solution  is  not  the  same 
as  pure  solvent,  and  (ii)  the  solvation  model  is  unsuitable. 
Fig.  3  shows  the  variation  of  change  of  net  charge  Ap(e)  as 
a  function  of  e  for  >C=0  and  CH3  of  PC.  As  e  increased, 


Table  l 


Equilibrium  constants,  Kd ,  Ks  and  critical  solvation  shifts  A  5' 
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Table  2 

Calculated  enthalpy  and  change  of  net  charge  in  vacuum  (e  =  1 )  and  in  dielectric  medium  (  eof  each  solvent) 
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Fig.  2  Illustration  of  solvation  shift  direction  and  change  of  net  charge,  Ap 
( in  vacuum  e  =  1 ),  for  PC:  (  ®  )  low  field  shift,  and  (  0  )  high  field  shift 
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Fig.  3.  Change  of  net  charge  A  p(e)  vs.  dielectric  constant  e  for  PC. 

A p(e)  for  > C=0  decreased  and  A p( e)  for  CH3  increased. 
At  e  around  3,  A  p(  e)  could  well  explain  A  8' .  But  for  linear 
carbonate  compounds,  decreasing  e  did  not  lead  to  a  good 
A p(e)  ratio  of  >C=0  to  CFL,  this  ratio  was  much  bigger 
than  its  A  8’  ratio.  For  cyclic  carbonate  compounds,  the 
>  C=0  oxygen  bonding  solvation  model  and  lowering  e  are 


suggested.  For  linear  carbonate  compounds,  another  solva¬ 
tion  model  is  needed  to  explain  the  A  8'  by  A  p.  By  consid¬ 
ering  these  factors,  A H(e)  may  show  better  correlation  with 
K s  ( it  must  be  related  to  A  G  =  —  RT  In  AT,  if  A  S  is  not  large ) , 
and  further  simulation  is  needed  in  this  sense. 


5.  Conclusions 

Characterization  of  a  pure  solvent  electrolyte  system  of 
cyclic  and  linear  carbonate  compounds  by  NMR  solvation 
shifts  and  MO  simulation  suggested  that  cyclic  compounds 
solvated  to  Li+  with  >  C=0  oxygen  and  the  dielectric  con¬ 
stants  of  their  electrolyte  may  be  smaller  than  those  of  pure 
solvents,  while  linear  carbonates  may  not  only  bond  to  Li  + 
with  >  C=0  oxygen  but  also  -O-  oxygen.  MO  simulation 
derived  a  net  charge  deviation  in  the  solvent  molecule  by  Li + 
solvation.  This  explained  the  higher  field  solvation  shift  of 
certain  atoms  in  solvent  molecules  even  in  the  case  of  bonding 
to  Li+  cation. 
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